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O ' Regular model sets, describing the point positions of ideal quasicrystallographic 

. tilings, are mathematical models of quasicrystals. An important result in mathemat- 

ical diffraction theory of regular model sets, which are defined on locally compact 
Abelian groups, is the pure pointedness of the diffraction spectrum. We derive an ex- 
tension of this result, valid for dense point sets in Euclidean space, which is motivated 
by the study of quasicrystallographic random tilings. 



1 Introduction 

An important question in mathematical diffraction theory concerns the problem which dis- 
tributions of matter diffract (Bombieri and Taylor (1986)). By now, there are only partial 
answers to this question which go beyond the crystallographic case. A mathematical ideal- 
ization of the set of atomic positions of a piece of matter are Delone sets (Lagarias (2000)). 
A subset A of M.'^ is called a Delone set if it is uniformly discrete and relatively dense. This 
means that there are radii r,R > such that each ball of radius r (resp. R) contains at 
most (resp. at least) one point of A. This class, however, is too general to obtain spe- 
cific results about spectral properties. It includes, for example, ordered structures such as 
crystals, as well as structures with disorder, and also amorphous systems. 

A special class of Delone sets are model sets (see for example Moody (1997)), which 
arise from a cut-and-project scheme. We will repeat their precise definition below. Model 
sets have strong regularity properties such as uniform point densities. They are mathe- 
matical abstractions of ideal quasicrystals, whose diffraction spectrum is experimentally 
known to consist of Bragg peaks only. By now, the study of model sets is a rather 
well developing subject, and their diffraction properties are well understood (Hof (1995), 
Hof (1998), Schlottmann (1998), Schlottmann (2000)). The most general description is 
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in terms of measures on locally compact Abelian groups, in generalization of Euclidean 
space (Sclilottmann (1998)). It is known that regular model sets are pure point diffrac- 
tive. A first proof of this fact was given by Hof (1998) and, in a more general setup, by 
Schlottmann (1998) , where they used a dynamical systems approach together with an 
argument due to Dworkin (1993). Recently, an alternative proof has been given by Baake 
and Moody (2002) , who considered Delone sets with certain additional properties, which 
include model sets. For these sets, they explicitly constructed a cut-and-project scheme 
and were able to prove the pure pointedness of the diffraction spectrum directly. As they 
showed, their results have a natural interpretation in terms of the theory of almost pe- 
riodic measures (Gil de Lamadrid and Argabright (1990)). An important assumption of 
the approach is the (uniform) discreteness of the point set A, arising from the physical 
motivation. 

In this paper, we will show that the diffraction result for regular model sets remains 
valid for certain dense point sets in Euclidean space, which appropriately generalize model 
sets. These point sets are not a mathematical curiosity but admit an interpretation in terms 
of averaged structures derived from random tilings with quasicrystallographic symmetries 
(Henley (1999), Richard et al. (1998), Richard (1999)): Soon after the discovery of qua- 
sicrystals, it became clear that there are two competing models for the description of their 
unusual diffraction properties, namely the ideal tiling model and the random tiling model. 
Whereas the ideal tiling model (i.e., a model set) leads to a diffraction spectrum consist- 
ing of Bragg peaks only, in the random tiling model a continuous background in addition 
to Bragg peaks was expected, as has been argued by Henley (1999). There is at present 
no rigorous treatment of diffraction properties of quasicrystallographic random tilings in 
d > 1, apart from the comparatively simple situation of disorder arising from independent 
random variables (Baake and Moody (1998)) and the investigation Kiilske (2001). For 
rigorous results about diffraction of crystallographic random tilings in d = 2 and random 
tilings in d = 1, see Baake and Hoffe (2000), Hoffe and Baake (2000). The derivation of 
Henley's prediction rests on non-rigorous arguments and uses an averaged point distribu- 
tion, which can be regarded as a generalized model set, being dense in Euclidean space. 
This distribution was assumed to exist in dimensions d > 2 (Henley (1999)). In the follow- 
ing, we give a rigorous exposition of a main part of this approach, which is, on the other 
hand, a generalization of the cut-and-project scheme for regular model sets to dense point 
sets (see also Hoffe (2001), which introduces the setup used in this paper, and the review 
article Baake et al. (2002), where the diffraction formula of Theorem ITUl is cited). After rec- 
ollecting the basic terminology and central theorems of diffraction theory of regular model 
sets in Section 2, we will define generalized model sets in Section 3 and show that they are 
pure point diffractive. This leads to a generalization of Poisson's summation formula to 
certain dense point sets. Section 4 will apply the results within Henley's framework for the 
description of diffraction properties of quasicrystallographic random tilings. A discussion 
of open questions and possible future work concludes the paper. 
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2 Diffraction of weighted regular model sets 



In this section, we recollect and extend results from diffraction theory of weighted regular 
model sets, which we shall generalism in the following section to the setup of a dense Dirac 
comb. 

Throughout the paper, we consider the situation of a cut-and-project scheme (Moody 
(1997)) with a Euchdean space E*^ and a Euchdean space R"*, called direct space and 
internal space: 

Let L C X M.^ be a lattice and denote the volume of its fundamental domain by 
|det(L)|. If TTi and 7r2 are the natural (orthogonal) projections of x onto R'' and 
M"*, respectively, then ni\iis assumed to be one-to-one and 7r2(-L) is assumed to be dense. 
Set L — 7ri(L) and let ()* : L — > 7r2(L) denote the mapping 7r2 o (ttiI^;)"^, also called the 
star-map. This is summarized in the following diagram. 



TTi 772 



u \ u y 

1-1 dense 
L L ^ {{x,x*)\x e L} 



(1) 



For an open and relatively compact subset W C W^, also called window, define the model 
set A{W) by 

A{W) := {x e L\x'' eW}. (2) 

Denote the volume of W by vol(W^). Model sets are Delone sets, i.e., they are both 
uniformly discrete and relatively dense (Moody (1997), Schlottmann (1997)). The model 
set A{W) is called regular ii $ W — int(W^) is compact and if dW has zero Lebesgue 
measure. 

Associated with the lattice L is its dual lattice (L)*, defined via 

(L)* = {y eM.'^ \ y X e Z ior all X e L}. (3) 

Denote its image in by L* = 7ri((L)*). A remarkable property of cut-and-projcct 
schemes is duality (Moody (1997), p. 418). We can dualize the given cut-and-project 
scheme to obtain a cut-and-project scheme for the dual lattice. By identifying direct and 
internal space with t^eir corresponding duals, we can see that the projection tti, restricted 
to the dual lattice (L)*, is one-to-one, and that the dual lattice has dense image 7r2((-L)*) 
in R"*. The corresponding star-map is defined on the set L*. We denote it by ()* again. 

Regular model sets have a well-defined density. Throughout the paper, we will consider 
infinite volume limits to be taken on sequences of balls. Existence of these limits may 
however be derived for more general van Hove sequences, see Schlottmann (2000) for the 
example of the density formula. Let S„(a) denote the closed ball of radius n centered at 
a e R'^. We set S„ = S„(0). 
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Theorem 1 (Density formula (Schlottmann (1998))). Let A{W) be a regular model set. 
Then 

,-»vol(S.(0)) l^.,^,^t;nB,,.) ) Idet(i)r 

uniformly in a and in u. □ 

This result may be used to consider sums over weighted regular model sets. Consider 
first the situation where the weight of a point x e A(W) only depends on x. 

Lemma 1. Let A{W) be a regular model set, / : M*^ — > C bounded and \x\'^^^~^'^\f{x)\ < C 
for some constants C > and a > 0. Then, for a and u e , the sums 

s{u,a):^ X] -^(^^ 

xeA(W+u)+a 

are absolutely convergent. Moreover, they are uniformly bounded in u and in a, the bound 
being proportional to vol(H^). 

Proof. For e N, consider the sums 

s<{u,a,R):= ^ |/(x)|, s>{u,a,R):= ^ |/(a;)|. 

xeh{W+u)+a x&K{W+u)+a 
\x\<R \x\>R 

In the density formula, the sequence on the Ihs is certainly bounded by twice its limit for 
almost all r. This implies for the number of points within a ball of radius n 

I (A(W^ + m) + a) n5„(0)| = \MW + u) r} Bn{-a)\ < cvo\{W)n'^ (n > no{W)), 

uniformly in u and in a, where c = 25*^/1 det(L)|, and Sd C M.'^ denotes the volume of the 
unit ball. The number uqCW) is increasing with decreasing volume of W. Together with 
/ bounded, the above estimate implies that s<(m, a, R) is uniformly bounded in u and in 
a with a number proportional to vol(l^). We derive a uniform bound on s>{u, a, R). By 
assumption, we have |x|'^+-'^'*""|/(x)| < C. Define gn ■— '^n<\x\<n+i 1/(^)1 estimate gn 

by 

< I (AiW + u) + a)n B^MO)\ < cvol(^)(n + 1)^-^ < ^^^^^ (4) 

for n > ni(W). Thus for R > ni(W) 

\s>iu,a,R)\ <J29n< 2cCyo\{W) ^ < oo, (5) 

n>R n=R 

since by assumption a > 0. The bound is independent of u and of a. Note that we have 
hmij_oo s> {u, a, R) ^ 0. □ 
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Remark. The above sums s{u, a) are absolutely convergent and bounded uniformly in u 
and in a under the milder assumption | < C for some constants C > and a > 0. 

This follows already from the uniform discreteness of A(W^) by a standard estimation. The 
additional property that there exists a bound proportional to vol{W), which we will use 
extensively below, is a result of the density formula. It would be interesting to consider 
whether the assumptions on / in the above lemma may be weakened, since the estimation 
in equation seems rather crude. 

Consider now the situation where the weight of a point x G A{W) depends only on 
its internal coordinate x*. This leads to Weyl's theorem on uniform distribution in the 
context of regular model sets (Kuipers and Niederreiter (1974), Schlottmann (1998)). 

Theorem 2 (Weyl's Theorem for regular model sets (Baake and Moody (2000))). 
Let A{W) be a regular model set, with compact, Riemann measurable W C M™. Let 
f : C be continuous with supp(/) C W . Then, for all a G Mf^, 

1 1 r 

™^°l(^^'(0)).eA(W(a) |det(L)|yH^ 

uniformly in a. □ 

The fundamental object in diffraction theory of weighted model sets is the weighted 
Dirac comb u defined by 

^= E /(^')'^- (6) 

x£A(W) 

where sup^^^^rJ\f{x*)\ < oo. This defines a complex regular Borel measure on M.'^, which 
is translation bounded, since A{W), being a model set, is uniformly discrete. Recall that a 
measure u is translation boundediS, for all compact K C W^, supyfz]^d\u!\{y + K) < Ck < oo 
for some constant Ck which only depends on K. Here, |ci;| denotes the total variation 
measure and y + K = {y + x\x G K}. In the following, we assume that the function 
/ satisfies the assumption of Weyl's theorem for regular model sets, i.e., we assume that 
/ : C is continuous with supp(/) C W. 

Diffraction properties can be expressed using the Fourier-Bohr coefficient cw{k) of the 
weighted Dirac comb uj, defined by 

Mk) = lim I Yl (7) 

r^oo vol i3r 

where A; G M*^ and a G R*^. We have the following theorem. 

Theorem 3 (Fourier-Bohr coefficients (Bernuau and Duneau (2000), Hof (1995))). 

Let K{yV^ be a regular model set, with compact, Riemann measurable W C M™. Let 
f : M™ — * C 6e continuous with supp(/) C W . Then, for all k G M.'^, the Fourier-Bohr 
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coefficient cw{k) exists and is independent of a e R*^. Its value is given as follows. For 
any k & L*, one has 

cw{k) = [ e'-^'*-f{u) du = -r^fl^i-kn, 

I det(L)| Jw I det(L)| 

and cw{k) — if k does not belong to the lA-module L*. □ 

Remcirk. In Bernuau and Duneau (2000), the theorem is proved only for the case where 
f{u) equals unity. The statement can be generalized to the situation described above by 
the same methods which lead to Weyl's formula, generalizing the density formula. For 
/c e L*, the theorem is a direct consequence of Weyl's density formula. 

For regular model sets, the weighted density of points p exists, because Weyl's theorem 
for regular model sets implies 

1 1 /■ 

p lim — , , ^ . uj{Br) = / f(u) du. (8) 

™vol(5,) ^ |det(L)| Jw 

This identity may be viewed as a particular normalization of admissible functions /, which 
we employ in the following. 

Diffraction is described by properties of the Fourier transform of the autocorrelation, 
which we now define (see also Hof (1995), Baake, Moody and Pleasants (2000)). Set 

u;{f) = uj{f), where f{x) — f{—x). Define truncated Dirac combs a;„ = ^\b„ 
cUn — i^^nT- The finite autocorrelation measures 

arc well defined, since c<j„ has compact support. Recall that the convolution of two measures 
p, V is defined as ^ * z/(/) = /jgd^E'* ^ dp{x) dv{y), being well-defined if at least one 
of them has compact support. The finite autocorrelation measures read explicitly 

x—y=z 

where A = A{W) — A{W) is the set of difference vectors of A (14^). It can be shown that 
the vague limit n — > oo leads to a unique autocorrelation 7^. 

Theorem 4 (Autocorrelation). Let A{W) be a regular model set, with compact, Rie- 
mann measurable W C . Let f : — > C 6e continuous with supp(/) C W . Then, the 
vague limit n ^ 00 of the finite autocorrelation measures 7^"''' leads to a unique measure 
■y^, called natural autocorrelation, being a translation bounded, positive definite pure point 
measure. It is given explicitly by 

7. = 5] v(z)5,, viz) = — -L- / f(u)f(u - z*) du. (11) 

I det(L)| Jwn{w+z*) 

□ 
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Recall that a measure /x is positive definite iS ^{g *g) > for all compactly supported 
continuous functions g. 

The proof of the theorem proceeds as follows (Baake, Moody and Pleasants (2000)). All 
autocorrelation coefficients ri{z) = \imn^ooVn{z) exist due to Weyl's theorem for regular 
model sets. They are locally summable, since A is closed and discrete. Then, defines 
a distribution over the space of all C°°-functions of compact support. The translation 
boundedness is inherited from u (see Hof (1995), Prop. 2.2). Finally, 7^ can be written as 
a certain volume-normalized convolution, which implies that it is a distribution of positive 
type. ^ 

Due to Bochner's theorem (Reed and Simon (1980), p. 331), the Fourier transform 7^^ 
of 7^; is a positive measure and also translation bounded (Baake, Moody and Pleasants 
(2000)). Recall that the Fourier transform /i of a tempered distribution /i is defined as 
Jiif) = fJ'i'p) for all Schwartz functions (f, where <f{y) = J e~^'^*^'^(/9(x) dx. It can be shown 
that the Fourier transform of the autocorrelation is a pure point measure (Schlottmann 
(2000), Baake and Moody (2002)). 

Theorem 5 (Diffraction formula (Bernuau and Duneau (2000))). LetA{W) be a regular 
model set, with compact, Riemann measurable W C M™. Let f : C be continuous 

with supp(/) C W . The Fourier transform 7^ of the autocorrelation measure 7^^ is a 
translation bounded, positive pure point measure. It is explicitly given by 



where L* is the projection of the dual lattice into direct space, and cw{k) are the Fourier- 



Recall that a measure /x is positive iff fi{g) > for all compactly supported continuous 
functions g >0. 

Remark. The pure pointedness of has been shown in Hof (1998), Schlottmann (1998), 
Baake and Moody (2002). The explicit formula for the discrete part of 7^ was proved 
in Hof (1995), but appeared earlier at different places in the physical literature. In the 
context of deformed model sets, which include regular model sets as a special case, the 
theorem appears in (Bernuau and Duneau (2000)). The diffraction formula is often stated 
for unweighted Dirac combs, but also holds for weighted Dirac combs, as is seen by an 
approximation of the weight function / by step functions, analogously to the proof of Weyl's 
theorem for regular model sets in Baake and Moody (2000) using the density formula. 

3 Diffraction of dense Dirac combs 

We now extend the above results to the situation of a dense Dirac comb. The corresponding 
proofs will rely on the above results. In the following, we assume a cut-and-project scheme 




Bohr coefficients of Theorem\^ 



□ 
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as in the previous section with the additional property that the canonical projections tti 
and 7^2, restricted to L, are both one-to-one, and the images of L are both dense in M.'^ 
and M"^, respectively. The star-map is then a bijection between L and L*. This additional 
assumption makes it possible to regard subsets B gM.'^ as windows, leading to model sets 
in internal space M*". 

Define the weighted Dirac comb uj by 



UJ = 

Since the Z-module L is dense in Euclidean space, the above sum is now well defined only 
under special assumptions on the weight function /. 

Theorem 6 (Weighted Dirac combs with dense support). Let f : M™ ^ C be 

bounded, and |a;|'^"^^"''"|/(a;)| < C for some constants C > and a > 0. Then, the 
weighted Dirac comb uj is a translation bounded measure. 

Proof. Let a compact i^' C M*^ be given. Cover K with a finite number of translated unit 
balls Wi such that K C ljr=i ^« ~ y ^ total variation measure of u is 

bounded by 

\u;\{K + y)< J2 l/(^*)l= E \f{x*)\=--aK{y). 

xGLn{W+y) x*€{Ln(W+y))* 

The additional assumptions on the cut-and-project setup, introduced at the beginning of 
this section, imply that (L fl {W + y))* is a regular model set with window (W + y) C M°'. 
Lemma H then yields that axiy) < bx < oo, uniformly in y, for some bx- This implies 
that is a measure, and that uj is translation bounded. □ 

We derive a generalization of Weyl's theorem on uniform distribution to dense point 
sets. 

Theorem 7 (Weyl's Theorem for dense point sets). Let f : M™ ^ C be continuous, 
and |a;|'^+^"''"|/(x)| < C for some constants C > and a > 0. Then 

1 1 r 
lim — , , ^ , > fix*) = / f(u)du, 

uniformly in a. 

Proof. For s G N, let Bs C M™' denote the ball of radius s centered at 0. The idea of the 
proof is to approximate the dense set L by the model sets A{Bs). Let Xs '■ I^™ [0, 1] be 
continuous with Xs{x) = 1 for |x| < s — 1 and Xs{x) = for |x| > s. Define the numbers 



Wr.. = V iXs ■ f){x*), W'^ = ^ [ iXs ■ f){u) du, 

^TT^r^ , ^ det L JB^ 



xeA{Bs)nBr{a) 
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We have linir 



lim^^oo w'g = w < oo. This follows from Lebesgue's dominated convergence theorem, since 
the functions \{xs ■ /)(^)| are bounded by uniformly in s, and f{u) is by assumption 

absolutely integrable (the integral being finite). 

We now show that lims^oo 'Wr,s = Wr uniformly in r and in a. We have the estimate 



w'g < oo due to Weyl's theorem for regular model sets. We also have 



\Wr — Wr s\ 



< 



vol (5^(0)) 



vol(S^(0)) 



vol (5^(0)) 



x&LnBr{a) 
( 



X 



x&A{Bs)r\Br{a) 



X 



\ 



x*&{Lr\Br{a))*\Bs 



x*e(LnB,.(a))* 
s-l<|a;*|<s 



1/(^^)1 



.X* (^(Lr\Br{a)Y\Bs^i 



We estimate the last term using equation (jSJ in Lemma [T] Since rii^Br^) < ni{Br^) for 
r2 > Ti, we get for s > ni{Bi) and r > 1 



\Wr — Wr,s\ ^ 



— ^ ( 2cCvol(5,) V -^1 . 
vol(5.) 1^ 



Since the bound is independent of r and a and vanishes as s ^ oo, the assertion follows. 

We now use a 3e argument to show that Wr — > w. Fix e > 0. We have — < e 
for s > sqIs). We also have \wr^s — < s for s > Si{e) uniformly in r. Lastly, we have 
\w'g — Wr^s\ < s for r > tq^e, s). Take s > max{so(e), Si{e)}. For r > tq^e, s), we thus have 



\w — Wr\ < \w — w'g\ + \w'g — Wr^s\ + I'U^r.s " | < 3 £. 

This establishes the convergence result of the theorem. 



□ 



For applications to diffraction, we restrict the class of admissible functions / in the 
remainder of this paper. Let / : M"^ ^ C be continuous, and |/(a;)| < C / \x\"^^'^^°' for 
some constants C > and a > 0. We can argue as in the previous section, using Weyl's 
theorem for dense point sets, that the weighted density of points exists. 



P ■-- 



lim 



1 



vol (5,. 



f{u) du. 



(13) 



det(L)| 

This identity may be viewed as a particular normalization of admissible functions /. Given 
a weighted Dirac comb uo with dense support, we consider its Fourier-Bohr coefficients. We 
define the finite volume approximations 

1 



Cr{k) 



vol(5^(0)) 



x&LnBr{a) 



X e 



-2'iTik-x 



(14) 
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where k and a G M''. The numbers Cr{k) exist (and may depend on a), as is seen by 
an argument similar to that used for proving the existence of the Dirac comb u. In the 
limit r — > CX3, we have the following result. 

Theorem 8 (Fourier-Bohr coefficients for dense point sets). Let f : ^ C be 

continuous, and < C for some constants C > and a > 0. Then, for all 

k G M.'^, the Fourier-Bohr coefficient c{k) = lim^^oo Cr-(fc) exists and is independent of a. 
Its value is given as follows. For any k & L* and for any a G M.'^, one has 

c{k) = [ e2-^"V(«) du = ^-^/(-r), 

|det(L)|yR- |det(L)| 

and c{k) = if k does not belong to the Z-module L* . 

Proof. The proof is analogous to the proof of Weyl's theorem for dense point sets above. 
For s G N, let C denote the ball of radius s, centered at 0, and define the numbers 




vol(5.(0)) ^ iXs-f){xne 



J^\lBSXs-fKu)e^-''-du, keL* 

0, k^L* ' 

0, k^L* ' 

We have lim^^oo Cr,s(fc) = c'^ik) < oo due to Theorem |21 We also have hm^^^oo c'^(^) = 
c'{k) < cxD by Lebesgue's dominated convergence theorem. We show that lim^^oo Cr,s(A;) = 
Cr{k) uniformly in r and in a. As in Theorem[7l we can then use a 3e argument to show that 
Cr{k) converges to c'{k) as r ^ oo. Thus c{k) = lim^^oo Cr{k) = c'{k), and the assertion of 
the theorem follows. 

For every G M'^, we have the estimate 



\Cr{k) - Cr,s{k)\ < 



x*&{LnBr(a))*\Bs x*e{LnBr{a))* . 

\ s-l<\x*\<s / 



vol (5^(0)) 



This is the same expression as in the proof of Theorem [7| We can repeat the previous 
argument and derive a uniform bound in r and in a, which vanishes as s ^ oo. □ 



The finite autocorrelation measures of the dense Dirac comb u, 



7i"^=^-n*^„, (15) 
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are well defined, since Un has compact support. They read explicitly 



VoUBn) 

zGL ^ x,y£LnB„ 

x—y=z 

We show that the limit n ^ oo leads to a unique autocorrelation. As a first step, we show 
that the pointwise limit exists. 

Lemma 2. Let rin{z) denote the finite autocorrelation coefficients of the weighted Dirac 
comb uj. The numbers rin{z) have a well-defined limit t]{z), which is a positive definite 
function. It is given by 

1 



lim 7]n{z) = r]{z) = ^ / f{u)f{u - z") du. (17) 
'^^'^ I det(L)| Jr™ 

Proof. This is an application of Weyl's theorem for dense point sets. We have 

x-zeLnBn 

= — TTTTT V fix*) fix* - Z*) y f(x*)f(x* - Z*). 

The first term in the last line converges to ri{z) by Weyl's theorem for dense point sets, 
and the second one, which we denote by r„(z), converges to zero, as we now show. Note 
that \x\ <n — \z\ implies |x — z| < n. Consider for n > \z\ the estimate 



1 



E \fix*)\\fix*-z*)\. 



V0\(Bn) 

^ :E*e(Ln(B„\B„_|,|))* 

This is a summation over a regular model set with window i?„ \ Bn-\z\- Due to Lemma ^ 
the sum is bounded by the volume of the window, which is a shell of thickness l-z] of the 
ball of radius n. Thus, the last expression vanishes like 1/n as n oo. The limit r]{z) is a 
positive definite function, since it is the limit of the positive definite functions rin{z). □ 

Theorem 9 (Autocorrelation). Let f : ^ C be continuous, and |a;|"*"'"^"'""|/(a;)| < C 
for some constants C > and a > 0. The weighted dense Dirac comb uj has the unique 
autocorrelation 

l^ = y^viz)5z, viz) = \~ / fiu)f{u- z*)du, 

det L JR'" 



zeL 



where 'j^j is a translation bounded, positive definite pure point measure. 
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Proof. Consider the regular model set autocorrelation measures 



7u,s = Y] Vs{z)Sz, Vsiz) = - / ixs ■ f)iu){xs ■ f){u - z*) du, 

det L JB.,n(B„+z*) 



where Bg denotes the ball of radius s G N, and = A{Bs) — A{Bs). Note that C 
As_|_i C L and set ris{z) = if z ^ A^. Consider the associated finite autocorrelation 
measures 



7/ 



z&As ^ ' x,y£A{Bs)nB„ 

x—y=z 



Since these measures arise from regular model sets, we have 7(i"] '-)^^s vaguely. The 
measures -^i^^s converge to 7^^ vaguely, as follows from Lebesgue's dominated convergence 
theorem: Note that ris{z) r]{z) and \ri,{z)\ < \ri\{z) := J^n. \f{u) \ \f{u - z'')\du. 

Note that 7|^| := J2z€L \v\{^)^z defines a measure: For 7|^| and compact K, we have the 



estimate 



zeLnw z*(^{Lr\W)* 



det(L)| 



/ I/Ml E l/(«-^')l du<^, 



,2*e{Lnvl/) 




where W = ljr=i ^« ^ covering of with a finite number of translated unit balls 
Wi. In the last equation, we exchanged summation and integration, which is justified 
by Lebesgue's dominated convergence theorem, since the terms in brackets are bounded 
uniformly in u due to Lemma ^ For g continuous with support on K, we then have 



lim 'y^Aa) = lim E ^«(^)^(^) = E ^™ Vs{z)g{z) = E v{z)g{z) = 7^(^). 

z&LnK z&LnK z&LnK 

We will show below that ^yj^l — > 7(1"^ vaguely, uniformly in n. The above results may then be 
used with a 3e argument as in Theorem [7| to show that 7!"^ 7^^. Since u is translation 
bounded, we further conclude that all finite volume approximations 7(i"^ are uniformly 
translation bounded (Hof (1995), Prop. 2.2), hence the limit is translation bounded. As 
it coincides with the pointwise limit, it is unique. The measures 'j!^'' are positive definite 
by construction. Since the positive definite measures are closed in the vague topology, 7^ 
is a positive definite measure. The explicit form of the vague limit shows that 7^^ is pure 
point. 

To show the uniform convergence of 7ii"], consider for a compact set K C 'M.'^ the 
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estimate 



vo\{B„) 



— 



E E iXs ■ f){x^){Xs ■ f){x* - z*) 



zeLnK xeLnB„ 

x-z£LnBn 



z£LnK x(^A{Bs)nB„ 

x-zeA{Bs)nBn 



< 



vo\{B„) 



vol(5„,) 



E E /(^^)/(^^ - - E (xs ■ mx^ixs ■ f){ 

zeLnK \^x<^LnBn xeA{Bs)nB„ 
( 



X — z 



zf^LnK 



\ 



E /(^^)/( 



X — 2; — 



E bCs-JW){Xs-f){x^-Z^) 



xeLnB„ 
\x-z(fLnB„ 



xeA{Bs)nB 

x-zfA(Bs)nBn 



J 



Both terms in the last inequahty may be estimated by 



1 



vol(5„) 



E 

z£LnK 



E 1/(^^)1 i/( 

x*e{LnB„)*\Bs 



X — z 



+ E 1/(^^)1 1/( 



X — z 



x*&(LnB„)* 
s-l<|3;*|<s 



\ 

I 

/ 



by the same argument as in the proof of Theorem [7| We may now interchange the sum- 
mations and use the fact that the sum over z is bounded uniformly in x* due to Lemma 
^ The remaining term is of the same form as that in Theorem [7| We conclude that it 
approaches zero as s — 00 uniformly in n. Thus 7^"] 7i"'' vaguely and uniformly in n, 
which concludes the proof. □ 

It now follows from the theorem of Bochner-Schwartz (Reed and Simon (1980), p. 331) 
that the Fourier transform of the autocorrelation is a positive, translation bounded mea- 
sure. An explicit expression for the discrete part (7^)pp of which is given in the following 
theorem, can be deduced from Hof (1995), Thm 3.4. The theorem states that the Fourier 
transform of the autocorrelation measure has no continuous component. 

Theorem 10 (Diffraction formula). Let f : ^ C be continuous and, for some 
constants C > and a > 0, |a;|''"+"'^^"|/(a:;) | < C. The Fourier transform %j of the 
autocorrelation 7^ of the weighted dense Dirac comb uj is a positive, translation bounded 
pure point measure and explicitly given by 



11. 



det(L)|2 



E i/(-fc*)r4, 



where the Z-module L* is the projection of the dual lattice. 
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Proof. We showed in the proof of Theorem El that the autocorrelation measures 7^,5 con- 
verge vaguely to 'j^^. Since the Fourier transform is continuous in the vague topology, the 
Fourier transforms of the autocorrelation measures 7^ converge vaguely to 7^. We will 
show that the vague limit leads to the above expression. 

Take a compact set K (ZM."^ and a covering W = [J^^i Wi of K with a finite number of 
translated unit balls Wi. Since (L* fl Wy is, by duality, a Delone set, we may order the 
numbers k E L* r\K to obtain a sequence {ki)i^^ with \k*- \ > \k*\ for j > i. Moreover, since 
Cg(fc) — ^ c{k) (in the notation of the proof of Theorem |H)), we may choose a subsequence 



such that \c'g.{ki)\^ — \c{ki)\ 



< 2~* on L* n for all j E N. For explicity of 



notation, we will suppress the index j in the following. Fix e > 0. We have 

00 00 00 00 



i=m+l 



i=m+l 



i=m+l 



i=m+l 



for m > mQ{e). The last sum can be made arbitrarily small, since it arises from the discrete 
part of the measure 7^, see Hof (1995), Thm3.4. Moreover, we have ^^^^ (|Cs(fci)P — |c(A;j)p) < 
e for s > SQ{e,m), since |c'^(A;)P |c(/;;)P uniformly in k. Fix now m > mo(£:). For 
s > so{e,m), we have 



1=1 



i=l 



< E \^'sm'+ 



i=m+l 



E(i4(fc.)p-ic(fc.)p) 



1=1 



E ic(^^)i'<2^. 



i=m+l 



For g continuous and compactly supported on K, this implies 

7:(^7) = lim 7:::;(^?) = hm e i4(^)I'^7(^)= E 

s— >oo s— >oo ' ' ' ' 

Since K was an arbitrary compact set, our claim follows. 



□ 



The above results lead to a generalized Poisson summation formula. Recall that, for a 
lattice L C M°', the Poisson summation formula is (Cordoba (1989)) 



det(L)| 



(1^ 



where | det(L) | is the volume of the fundamental domain of the lattice L, and the sum ranges 
over all points of the dual lattice L*. For regular model sets A(iy) C M'^, a generalized 
Poisson formula is given by 



(19) 



as is readily inferred from Theorem El This identity is not symmetric because the first sum 
ranges over the uniformly discrete set A = A(l^) —A{W), whereas the second sum ranges 



14 



over the projection of the dual lattice L*, which is a dense set. For weighted Dirac combs 
defined on a dense Z-module L C M'', we obtain from Theorem^] the formula 



where the sums on both sides are defined on a dense set. Within this symmetric setup, 
it may be easier to investigate which constraints on an underlying point set an identity of 
the above type imposes (Lagarias (2000), Problem 4.1). 

4 Applications to random tilings 

In this section, we explain how the above results can be used to infer diffraction properties 
of random tiling ensembles. 

A tiling of M'^ is a face-to-face space filling with tiles from a finite set of M'^-polytopes 
called prototiles, without any gaps and overlaps (Richard et al. (1998)). There might be a 
number of additional packing rules specifying the allowed configurations. Associated with 
a tiling is the set A of all vertex positions of the tiling. Take a ball Bn C M'^ of radius 
n and count the number of different patches of Bn, where we identify patches which are 
equal up to a translation. If the number of allowed patches increases exponentially with 
the volume of the ball, we call the set of all tilings a random tiling ensemble. 

The usual ideal quasiperiodic tilings like the Penrose tiling, the Ammann-Beenker tiling, 
and others, have the special property that the vertex positions of its tiles form a (regular) 
model set. Essentially, the set of all tilings is given by the collection of the original model 
set together with model sets of arbitrarily translated window. For tilings of Mf^ whose vertex 
positions constitute a Delone set and are described by a primitive substitution, it can be 
shown that the number of allowed patches of radius n grows asymptotically proportional 
to n'^ (Lenz (2002)). 

Relaxation of the packing rules for these tilings usually results in a random tiling 
ensemble with strictly positive entropy (Henley (1999)). Since entropy is an indication 
of disorder, one should expect a non-vanishing continuous component in the diffraction 
spectrum in addition to a discrete part. By construction, the set of vertex positions of 
each random tiling can be mapped into internal space via the star-map, resulting in a 
distribution, which may be supported on an unbounded domain in internal space. This is 
different from the distribution of a model set A{W), which is the characteristic function 
Iw of the window. 

A natural object to consider is the averaged distribution, where we take the average 
over all random tilings. Here, we adopt the normalization that G A for all random tilings 
A. The averaged distribution defines a weighted point set of M.'^ which is supported on a 
countable, dense subset of M'^. 

Henley (1999) gives non-rigorous arguments that the averaged distribution (for infinite 
tilings) exists in dimensions d > 2. In dimensions d < 2, the width of the averaged 




(20) 
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distribution for finite patches on balls should diverge with the radius n of the ball. In 
d = 1, the divergence is with the square root of the system size, in d = 2, the divergence is 
logarithmic. He concludes that in dimension d < 2 the diffraction spectrum of a random 
tiling should have a trivial discrete part. 

In d = 1, his arguments can be made rigorous (Baake et al. (2002)). Consider, for 
example, the Fibonacci model set. Other tilings may be treated similarly. Let Z[r] = 
Z + Zr C M, where r = (1 + V5)/2 is the golden mean. The ring Z[r] is the ring of integers 
of the quadratic field Q(t) = Q(v^). Let ()* denote the automorphism of Q{t) that maps 
y/b I— »• —Vb. The set L = {(x, x*) | x G Z[r]} is a lattice in M x M. The Fibonacci model set 
is the set A(W) = {x E Z[r] | x* G (—1, r — 1]}. Consecutive points of the model set have 
distances 1 or r. If we regard the point positions as left endpoints of half-open intervals 
u, V of length 1 and r, we get the Fibonacci tiling. A Fibonacci random tiling is an arbitrary 
sequence of intervals u,v such that the frequency of m is a.s. equal to 1/r. The averaged 
distribution for a patch of size is, to leading order in A^, given by (Hoffe (2001), Baake 
et al. (2002)) 

where erfc(x) = ^ e~*^ dt denotes the complementary error function. The distribution 

width therefore grows with the system size as a/TV. Together with the above analysis, this 
may lead to a trivial Bragg peak at the origin a.s N —>■ oo. This corresponds to the behavior 
of the Fibonacci random tiling, whose diffraction spectrum has been computed explicitly 
(Baake and Hoffe (2000)). It consists of an absolutely continuous component in addition 
to a trivial Bragg peak at the origin. It would be interesting to compute the continuous 
component of the averaged structure, which may be given by large corrections to the 
asymptotic behavior in equation ()21|). 

For d > 1, there are no rigorous results about the averaged distribution of quasicrys- 
talline random tilings, which is due to more restrictive matching rules for the prototiles. 
(In d = 1, the only matching rule is the face-to-face condition, resulting in Bernoulli ensem- 
bles, which are easy to analyze.) There are however numerical investigations for a number 
of quasicrystallographic random tilings in d = 2 and d = 3. 

In d = 2, the averaged distribution of the Ammann-Beenker random tiling appears to 
be of Gaussian type, with a distribution width diverging logarithmically with the system 
size (Hoffe (2001)). Thus, the situation is similar to the d = 1 case discussed above. A 
numerical analysis of the diffraction measure indicates a trivial Bragg peak at the origin, 
together with a singular continuous component (Hoffe (2001)). This behaviour is believed 
to be generic in d = 2 (Henley (1999)). 

Ind >3, the averaged distribution is predicted to exist (Henley (1999)). Our result then 
implies that the averaged structure is pure point diffractive, if the distribution function 
is well-behaved. However, the diffraction picture of a typical random tiling in c? > 3 
is expected to display an diffuse background in addition to Bragg peaks (Henley (1999)). 
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The averaging on the level of the Dirac comb thus has the effect of extinguishing continuous 
components of the diffraction spectrum. 

To conclude, it is necessary to prove Henley's statements about the existence of an 
averaged distribution with finite width in d > 3. Furthermore, it would be interesting to 
investigate the relation between diffraction properties of the averaged structure and the 
random tiling, in particular whether the discrete parts of both structures coincide and 
whether the continuous parts coincide in < 2. Concerning the results of this paper, a 
natural question is how they can be extended beyond the Euclidean case towards the setup 
of locally compact Abelian groups. (The proofs of the pure pointedness of the diffraction 
spectrum used the fact that is isomorphic to its dual.) Another aspect concerns the 
connection to the theory of almost periodic measures (Gil de Lamadrid and Argabright 
(1990)), which may also be used for a characterization of diffraction spectra with continuous 
components. 
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